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Abstract 

Thyroid hormone (TH) receptors (TRs) play central roles in metabolism and are major targets for pharmaceutical 
intervention. Presently, however, there is limited information about genome wide localizations of TR binding sites. Thus, 
complexities of TR genomic distribution and links between TRfi binding events and gene regulation are not fully 
appreciated. Here, we employ a BioChIP approach to capture TR genome-wide binding events in a liver cell line (HepG2). 
Like other NRs, TRfi appears widely distributed throughout the genome. Nevertheless, there is striking enrichment of TR(3 
binding sites immediately 5' and 3' of transcribed genes and TRfi can be detected near 50% of T3 induced genes. In 
contrast, no significant enrichment of TRfi is seen at negatively regulated genes or genes that respond to unliganded TRs in 
this system. Canonical TRE half-sites are present in more than 90% of TRfi peaks and classical TREs are also greatly enriched, 
but individual TRE organization appears highly variable with diverse half-site orientation and spacing. There is also 
significant enrichment of binding sites for TR associated transcription factors, including AP-1 and CTCF, near TR peaks. We 
conclude that T3-dependent gene induction commonly involves proximal TRfi binding events but that far-distant binding 
events are needed for T3 induction of some genes and that distinct, indirect, mechanisms are often at play in negative 
regulation and unliganded TR actions. Better understanding of genomic context of TR binding sites will help us determine 
why TR regulates genes in different ways and determine possibilities for selective modulation of TR action. 
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Introduction 

Previous unbiased studies of the genomic distributions of 
nuclear hormone receptors (NRs) have revealed several unexpect- 
ed aspects of their mechanisms of action [1-9] Early models 
suggested that NRs bind DNA response elements (REs) that are 
commonly located within target gene promoters. From these 
positions, NRs modulate nearby gene expression via hormone- 
dependent changes in coregulator binding that, in turn, lead to 
changes in chromatin organization, histone modification and 
RNA polymerase II recruitment and processivity [10,11]. ChIP 
studies, however, have revealed that: i) cognate NR REs are widely 
distributed throughout the genome and often far from obvious 
target genes, ii) ligand-dependent NR/DNA binding activities that 
were not seen in vitro can emerge in vivo and iii) NRs can display 
unexpected DNA binding preferences. Additionally genome-wide 
studies have confirmed suggestions that NRs can sometimes be 
recruited to target genes via interactions with heterologous 
transcription factors (TFs) and also revealed that NRs cooperate 
with unique subsets of heterologous TFs in composite modules, 
recognized through proximity of their DNA binding sites 
[5,7,12,13]. 

Currently, there are no comprehensive genome wide studies of 
thyroid hormone (TH) receptor (TR) binding sites in liver cells and 



few in other cell types. It is therefore likely that complexities of TR 
genomic distribution are not fully appreciated. TH acts through 
two closely related receptors (TRot and TRfi) which are both 
members of the NR family, to regulate genes involved in cell 
growth and differentiation, homeostasis and energy metabolism 
[14,15]. Both TRs commonly bind to specific TH response 
elements (TREs) within target gene 5' proximal [16]. Nevertheless, 
there are many examples in which TREs are located within, or 
downstream of, their T3 regulated target genes [17,18]. Further, a 
ChlP-on-chip analysis that interrogated TRfi binding events in the 
mouse cerebellum at locations from —8 KB/+2 Kb of known 
transcriptional start sites (TSS) revealed that 40% of 91 detectable 
TR binding events were in introns and not the proximal promoter 
[7] . Finally, a recent study, of similar design to one the described 
here, analyzed genomic localization of exogenously expressed TRs 
in a mouse neural cell line and showed that many genes that are 
induced by T3 lack a proximal TR binding site [19]. Thus, it 
seems likely that TREs are not restricted to proximal promoters 
and unbiased analysis of TR binding events is needed to define the 
extent of TRE locations and relationships to T3 regulated genes. 

Although predominant models of TR action suggest that they 
act similarly to other NRs, subtle differences in the nature of gene 
regulation by TR can only be addressed by a non-biased, genome- 
wide analysis. It is thought that TRs interact constitutively with 
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DNA as heterodimers with retinoid X receptors (RXRs) and this 
lack of ligand dependency has been confirmed for many individual 
[20,21]. There are, however, suggestions that TRs are released 
from subsets of response elements upon hormone binding [22], 
based on studies of TR DNA binding preferences in vitro and 
investigation of TR/DNA interactions within promoters of 
negatively regulated genes. Neither idea has been tested at the 
genome level. 

Although many TREs have been characterized, the defining 
components of TRE sequences have never been fully resolved. 
Canonical TREs that have been identified by conventional 
approaches are composed of single or multiple degenerate direct 
repeats of the AGGTCA half site spaced by four base pairs (DR-4) 
and comparison of these elements allowed establishment of a TRE 
consensus [23-27]. Analysis of individual TRE sequences has 
nevertheless also revealed significant individual variations in half 
site sequence and spacing and atypical half-site orientations, 
including everted and inverted palindromic repeats (ERs and IRs) 
[28] and there are suggestions that half site-proximal sequences 
may influence affinity for TRs versus other NRs that recognize 
AGGTCA half sites [29-32]. Further, TR-mediated transrepres- 
sion of thyroid stimulating hormone, thyrotropin releasing 
hormone and superoxide dismutase involve distinct elements that 
do not resemble [33-36]. 

TRs are known to interact with several heterologous TFs, but 
the extent of TR-TF interactions and identities of these proteins 
are poorly understood. TRs bind AP-1 [37,38], CREB [39] and 
p53 [40,41] and also form composite elements with the chromatin 
boundary defining protein CCCTC binding factor (CTCF) [42]; 
recent computational analysis suggested up to 18% occurrence of 
TREs proximal to CTCF binding sites [43]. While these TFs and 
others may play important roles in TR function, only a whole- 
genome binding analysis can reveal the extent of TR/TF binding 
site overlap and functional significance of TR/TF interactions. 

To begin to understand links between genomic organization of 
TR binding sites and TR dependent gene regulation, we 
characterized genome-wide TR binding events using a BioChIP 
approach [44]. Our findings point to different requirements for 
TRfS at T3 induced target genes versus repressed genes or genes 
that respond predominantly to unliganded TRs. We also observe 
unexpected features of TR binding site distribution around T3- 
induced genes, complex TRE composition and obtain evidence for 
TR cooperation with several heterologous TFs. Better under- 
standing of influences of genomic context upon TR action will 
help us determine why TR regulates different genes in different 
ways. 

Materials and Methods 

Reagents 

Sources of reagents are as follows: primers from Integrated 
DNA Technologies (Coralville, IA), TRfil antibody from Pierce 
Biotechnologies (MAI -2 16, Rockford, IL), acetylated histone H3 
K9 antibody from Cell Signaling Technology (#9671, Danvers, 
MA), M2-Streptavidin beads from Life Technologies (Carlsbad, 
CA), T3 from Sigma-Aldrich (St. Louis, MO), HepG2 cells from 
the American Type Culture Collection (Manassas, VA), Fugene 
HD from Roche Applied Science (Penzberg, Germany). 

Plasmids 

Reporter constructs were generated by mutagenesis of pGL4. 13 
[MinP] (Promega, Madison, WI), using primers containing TREs 
described in Table S 1 . The adm proximal promoter reporter was 
purchased from Switchgear Genomics (Menlo Park, CA). BIRA 



and BLRP-TEV vectors were provided by David Moore, Baylor 
College of Medicine, Houston, Texas. To generate BLRP-TEV- 
TRP, human TRpi cDNA (NM_000461) was PCR amplified 
with primers bearing sequences complementary to the multiple 
cloning region of BLRP-TEV vector (5'-aGCCTCGACGGTA- 
CCGATATCCTCGAgTGACTCCCAACAGTATGACAGAAA- 
ATGGCCTTA-3 ') [44,45]. Products were sequenced and insert- 
ed into the BLRP-TEV vector using a modified version of the 
standard Quickchange mutagenesis protocol (Stratagene, Santa 
Clara, CA). 

Cell Culture and Transfection 

TRP-BioChIP cells were created by stable transfection of TRfi 
and BiRA expression vectors and subjected to dual selection. 
HepG2 cells and derivatives were cultured in DMEM, with 
Penicillin and Streptomycin and subcultured every 2-3 days. Prior 
to transfection assays, cells were supplemented with resin-stripped 
FBS to remove T3 and transfected with constructs indicated and 
luciferase reporter activity was normalized by cotransfection of 
pRL which contains the open reading frame of renilla luciferase. 

BioChIP and ChlP-PCR 

BioChIP was performed as described [44]. Briefly, cells were 
plated in DMEM containing resin-stripped FBS, treated with 
1 00 nM T3 and crosslinked with 1 % Formaldehyde and 
chromatin fragmented with a Bioruptor sonicator (Diagenode). 
Binding and washing steps were performed as previously described 
[44] and resulting samples analyzed by QPCR with indicated 
primers (Table S2) or by sequencing using the methods indicated 
below. Samples were used in ChIP analyses were analyzed in 
triplicate, and three pooled replicate samples were used for ChlP- 
Seq analysis. 

Sequencing 

Solexa libraries were prepared using the NEBNext DNA 
Sample Preparation Kit (NEB) and Illumina PE adapters 
(Illumina) from three pooled replicates of control cell line and 
TRP-expressing cells +/— T3. Libraries were sequenced on a 
Solexa GAIIx following standard procedures. Solexa sequencing 
data was aligned to HG19 with Bowtie 0.12.7 [46]. Binding peaks 
were analyzed using QuEST 2.4 [47]. 

Sequence Data Analysis 

Genome-wide binding data was analyzed with utilities in the 
Cistrome Portal [48]. Patterns of genome-wide binding were 
analyzed by [48]. Sequence motif analysis was completed with 
MEME analysis [49] and custom Scripts (Table S3). Gene 
ontological analysis was completed with Gene Codis [50] . 

RNA Expression Analysis 

For arrays, RNA samples, including three replicates for each 
treatment group, were extracted using the Qiagen RNeasy kit 
(Valencia, CA). cRNA was prepared with the Totalprep96 kit (Life 
Technologies, Grand Island NY) and hybridized to Human HT- 
12 V4.0 gene expression beadchips (Illumina), according to 
manufacturers' instructions. Bead Intensity values were extracted 
by Genome Studio, and processed with the Lumi package of 
Bioconductor [51], by quantile normalization of log2-transformed 
expression values, with an offset value of 50; fold changes were 
calculated for each target, including correction for multiple 
hypotheses. Realtime PCR was used to analyze gene expression 
levels using the primer sets (Table S4) and standard SYBR Green 
reagents (Roche) on a Roche 480 Instrument. 
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EMSA 

Human TR(3 and RXRa DBD-LBDs were expressed as 
described below. Double-stranded oligonucleotides were labeled 
by standard polynucleotide kinase methodology and incubated 
with cellular extracts and 1 ug of poly(dl-dC) (Sigma), in a binding 
buffer containing 25 mil HEPES, 50 niM KC1, 1 mM dithio- 
threitol, 10M ZnS04, 0.1% Nonidet P-40, 5% glycerol. After 
30 minute incubation at room temperature, complexes were 
resolved on a 5% non-denaturing gel and visualized by 
autoradiography. For some experiments, we used infrared Dye 
EMSA in which TREs possess a 5 ' amine moiety on their forward 
strand and were labeled with 2X-excess of IRDye800-NHS-esther 
(Li-Cor Biosciences, Lincoln, NE) and purified by silica column. 
Complexes were resolved on a 5% non-denaturing gel and 
visualized on a Li-Core Odyssey Imaging System. 

Protein Purification 

RXRa or TRP DBD-LBD was expressed in a BL21 DE3 
Escherichia coli strain (Invitrogen), using plasmids bearing these 
ORFs downstream of a 6X his tag, under control of a T7 
promoter, based on the pet28a vector (Merck, Whitehouse 
Station, NJ). Protein expression was induced overnight at 18 
degrees Celsius with 0.2 mM Isopropyl-P-D-thiogalactopyrano- 
side (IPTG). Bacteria were pelleted and lysed by sonication in lysis 
buffer (20 mM Tris pH 8.0, 0.5 mM NaCl, 100 uM phenyl- 
methylsulfonyl fluoride (PMSF)). Protein was purified with a 
nickel-IDA agarose (Sigma-Aldrich St. Louis, MO). TR(3 was 
expressed in vitro using the T7 Quick-Coupled Transcription 
Translation System, according to manufacturer's instructions 
(Promega, Fitchburg, Wisconsin). 

Mice 

Indicated mouse liver samples were harvested in T-PER tissue 
lysis buffer with added protease inhibitor cocktail (Pierce 
Bioreagents) from 4 week-old C57 mice maintained on standard 
chow diet, sacrificed after overnight fast. All procedures were 
conducted with full approval of the TMHRI Institutional Animal 
Care and Use Committee. 

Results 

BioChlP-TR System 

We used BioChIP [44] to capture TR genome-wide binding 
events. We employed a dual-stable selection protocol to express 
human TRpi bearing an in frame N-terminal BLRP tag and the 
E. coli BirA biotin ligase in HepG2 hepatocyte cells (referred to 
hereafter as B7B cells). We isolated pools of transfected HepG2 
cells and verified TRP expression by real-time PCR (Fig. SI A). 
Selected cells displayed approximately eight-fold elevation of TR(3 
transcripts relative to low levels of TR(3 mRNA observed in 
control HepG2 cells that were stably transfected with BirA and 
empty vector (shown) and parental HepG2 cells (not shown). TR(3 
transcript levels were unaffected by T3 and there were no obvious 
changes in TRa transcripts. Western analysis confirmed that B7B 
cells expressed increased TRP protein versus control HepG2 cells 
(Fig. SIB). However, it was notable that TRP protein expression 
levels were lower than seen in equivalent amounts of mouse liver 
extracts indicating that this expression strategy did not result in 
supraphysiologic TRP protein levels. 

Exogenously expressed Biotin-tagged TRP displayed normal 
DNA binding and transactivation function (Fig. S1C-F). Gel shift 
analysis revealed increased TRE binding activity in extracts of 
B7B cells versus control HepG2 cells and that this complex 
supershifted with antibodies against TRP and RXRa (Fig. SIC, 



left panel), with an extent of superhshifting similar to that seen in 
previous studies [52,53], co-migrated with in vitro translated 
RXRa-TRP heterodimer and not complexes formed by in vitro 
translated TRP or RXRa alone (center panel) and displayed a 
modest mobility shift in the presence of T3 (right panel). Thus, 
TRP adopts standard heterodimeric form in these cells. Trans- 
fection of TRE-dependent luciferase reporters (DR-4 and ER-6) 
into the B7B cells revealed amplification of T3-dependent 
luciferase activity versus control HepG2 and this effect displayed 
T3 dose requirements that resemble prior results obtained with 
stably expressed TRP in this cell type [54] (Fig. SID). We also 
observed enhanced T3 activation of known direct human TR 
target genes, low density lipoprotein receptor (LDLR) [55], B-cell 
lymphoma 3 (BCL3 [56]) (Fig. S1E) relative to weak T3 induction 
seen in HepG2 cells that express BiRA alone (Fig. S1E) and 
parental HepG2 cells (not shown). For LDLR, this T3-dependent 
increase was also reflected in increased protein levels as judged by 
western analysis (Fig. S1F). Thus, overexpressed TRP displays 
normal heterodimer formation and transcriptional activity in this 
cell background. 

There was close concordance between patterns of T3 response 
in B7B cells and HepG2 cells that stably express Flag-tagged TRP 
and were previously created by our group [49] . Both the new B7B 
cells and Flag-TRP cells expressed similar levels of TRP transcripts 
(Fig. S2A) and exhibited very similar levels of T3 induction of 
endogenous target genes (Fig. S2B and not shown) and a standard 
DR-4 driven luciferase reporter (Fig. S2C). Thus, properties of 
B7B cells resemble those of other HepG2 cells that stably express 
TRp. 

Enrichment of TRB near transcribed genes 

We performed ChlPseq analysis of TRP DNA binding. Briefly, 
we treated B7B cells +/-T3, precipitated biotinylated BLRP-TRP 
and associated DNA fragments using a streptavidin-based 
purification system and analyzed TRP associated DNA sequences 
(Methods). Data were mapped to the human genome (hgl9) using 
Bowtie 0.12.7 [46], and genomic regions that were over- 
represented in TR-enriched samples, relative to unenriched 
control DNA, were assessed using the QuEST package, version 
2.4 [49]. 

Genomic distributions of TRs were similar to those of other 
NRs. We detected 5,791 TRP binding events in untreated TRP 
BioChIP cells and 6,792 in equivalent samples treated with T3 
(Fig. 1, S3 and Table S5); within the range of NR binding site 
representation in other ChlPseq studies. Distributions of TR DNA 
binding peaks were also similar to other NRs; around 15% of TRP 
binding events mapped to proximal promoters (within 5 KB 
upstream of transcriptional start sites, TSS) and large fractions of 
TRP peaks were located within introns or intergenic regions with 
the remainder associated with downstream region, 5 ' untranslated 
region (UTR), 3'UTR and exons (Fig. 1A). There were no 
changes in overall distribution of TRP +/— T3. 

Even though TRP is widely distributed throughout the genome, 
we observed striking enrichment of TRP near transcribed genes. 
We compared numbers of TRP binding events assigned to 
different genomic regions (Fig. 1A) versus percentage representa- 
tion of each region within total genome sequences (Fig. IB). This 
revealed a 1 2-fold enrichment of TRP binding within proximal 
promoters, 2-fold enrichment within 5 'UTR regions and 11 -fold 
enrichment at immediate downstream regions. Conversely, 
intergenic binding events were less frequent; around 50% of the 
genome is classified as intergenic and only around 40% of binding 
events mapped to these regions (Fig. IB). There was no 
enrichment of TRP binding events within exons or introns relative 
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Figure 1. Characterization of genomic binding events. A. Distributions of TRfj binding peaks across specific genomic regions in the absence 
(black) and presence (grey) of T3. B. Bar graph representing relative enrichment of TR(3-bound regions within genomic intervals specified. Gene- 
proximal regions, including promoter regions, 5'UTR regions and downstream regions were highly enriched in TR|3-bound regions of the genome. C. 
Frequency distribution plot of binding events in regions proximal to transcriptional start sites (TSS) +/— T3 (blue and red, respectively). The x-axis 
represents nucleotides upstream and downstream of the TSS, y-axis represents numbers of binding events. 
doi:1 0.1 371 /journal.pone.0081 1 86.g001 



to their overall genomic representation. The preference of TR|3 
for proximal promoters of transcribed genes was also reflected in a 
frequency distribution graph, which revealed significant enrich- 
ment of TR(5 binding around the TSS (Fig. 1C). Similar to overall 
binding site distribution, T3 did not affect enrichment of TRfi 
binding events near actively transcribed genes (Fig. IB) or the TSS 
(Fig. 1C). Thus, TR(5 binding events are widely distributed 



throughout the genome but enriched near transcribed genes. 
Moreover, T3 does not trigger large scale redistributions of TRf} 
between different functional regions of the genome. 

TR Binding Events are Associated with T3-lnduced Genes 

To understand relationships between TR binding events and 
target genes, we compared TRR binding with TRpVT3 dependent 
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changes in gene expression. We treated B7B cells with T3 for 
8 hrs; which allows induction of a large number of T3 induced 
genes through direct regulatory events in HepG2 cell backgrounds 
[49] . W e performed an array-based analysis of gene expression 
and detected 411 T3-regulated genes with at least 1.7 fold 
regulation (adjusted P<0.05), including 282 induced genes and 
129 repressed genes. Large numbers of genes also displayed 
changes in basal activation or repression in the presence of 
unliganded TR(5 versus control cells (TRP effect; Fig. 2A). 
Numbers and identities of T3- and TRP-dependent genes were 
consistent with previous studies of HepG2 [54,57]. 

There was a high prevalence of TR binding events near T3- 
induced genes. More than 10% of T3 target genes contained a 
TRP binding site within 1 KB of the TSS relative to <1% of 
unregulated genes. The percentage of T3 target genes associated 
with binding events increased as larger distances were considered; 
close to 50% of T3-induced genes displayed TRP binding within 
25 KB of the TSS whereas less than 10% of unregulated genes 
displayed binding events (Fig. 2B). Searches of the human genome 
for consensus TREs (DR-4, IP-6) using custom in house algorithms 
(supplementary material) found no evidence for enrichment of any 
of these elements near T3 induced genes versus their overall 
representation within the genome (not shown). Thus, TRP is 
enriched near T3 induced genes but this does not reflect clustering 
of classical TREs at these locations. 

By contrast, other classes of TR-regulated genes were not 
clearly associated with TRP binding events (Fig. 2B). We observed 
possible weak enrichment of TRP binding sites near genes that 
were activated by unliganded TRs (upper panel, TRP effect); 
around 15% were associated with nearby TRP binding sites versus 
10% of unregulated genes. More surprisingly, we did not detect 
obvious TRP enrichment near genes that were transcriptionally 
repressed in response to T3 or by unliganded TRs (Fig. 2B, lower 
panel; Discussion). 

Clusters of TR Binding Events near T3 Induced Genes 

We performed a detailed survey of TRP binding near genes 
regulated more than 2.5-fold by T3 in the microarray (Fig. 3 and 
supplementary data). As expected, large numbers of TRP binding 
events occurred near highly induced targets; >50% of these genes 
exhibit TRP binding within 25 KB and >10% within 1 KB. 
Genes in this TRP bound, T3 induced, set were physiologically 
relevant and representative of TRB's known functions in liver (not 
shown). There was high correlation between identities of highly T3 
responsive genes in B7B cells and observed T3 responses that were 
insensitive to cycloheximide (CHX) treatment in related HepG2 
cells that express flag tagged TRP (see data in reference [49] and 
not shown). Further, these genes displayed rapid onset of T3 
response in Flag-TRP cells at early time points (see [54]). This 
implies that this gene set is highly enriched for direct TR targets. 
By contrast, only one T3 repressed gene (PDE2A) displayed 
nearby binding events, within 1 Kb of the transcriptional start site, 
and no other repressed genes displayed detectable TRP binding 
events. 

Closer analysis of patterns of TRP binding near selected genes 
revealed clusters of TRP binding sites in diverse distributions. In 
Figure 4, TRP binding events are represented by blue bars 
defining the extent of sequences precipitated by the BioChIP 
experiment and red bars representing TRP peak positions defined 
by QuEST. We observed TRP binding near both verified human 
TRP target genes LDLR and BCL3 (Fig. SIC) and this mapped to 
the 5' region of the transcription unit for LDL-R and to the 5' and 
3' regions for BCL3 (Fig. 4A). Similar "5' only" and 5'+3' 
distributions were also seen for other genes (not shown and see Fig. 
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Figure 2. Links between TRp Binding and regulatory events. A. 

Bar graph representing numbers of genes that display positive 
regulation (upper panel) or negative regulation (lower panel) that 
met statistical significance and an arbitrary +/- 1.7-fold cut-off in an 
array-based analysis of TRB-BioChIP cells +/-T3 or in TRfS-BioChlP cells 
versus parental cells that lack TRp (THRB effect). B. Bar graph 
representing percentages of TRB binding events within 1 KB, 5 KB or 
25 KB of the TSS of T3 induced, TRB induced or unaffected genes 
(upper panel) or T3 or TRB repressed genes (lower panel). Progressively 
lighter shading in the bar graph columns represents increasing distance 
from the TSS. 

doi:1 0.1 371 /journal.pone.0081 1 86.g002 

S4). Interestingly, for both genes, TRP peaks detected by BioChIP 
were close to known TREs in the proximal promoter region 
[55,56]). We also observed other types of TRP binding site 
distributions near highly induced genes; these included intronic 
(NCOR2 and ADSSL1) and 3' only (SOX7). We confirmed TRP 
binding at predicted sites by conventional ChIP with an antibody 
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Figure 3. Patterns of TR|5 binding and transcriptional regula- 
tion. Heatmap depicting log2-transformed expression levels (left) and 
TR(3 binding events within 1 KB, 5 KB or 25 KB of TSS (right) of genes 
that met statistical significance and an arbitrary +/-2.55-fold cut-off of 
gene induction in TRp-BioChIP cells +/— T3. Columns reflect the average 
of three experimental samples. Expression values in heatmap are as 
indicated by color scale (bottom, green indicating —5.7-fold repression, 
red indicating 13-fold induction), and location of binding events within 
the indicated ranges are depicted by the presence or absence of black 
bars in the three right-most columns. 
doi:1 0.1 371/journal.pone.0081 1 86.g003 

specific for human TR|J (Fig. 4B, amplified regions are represent- 
ed by black bars in Fig. 4A) and verified that there was increased 
T3-dependent gene induction versus control HepG2 cells with 
qPCR (Fig. 4C). Thus, TR[J binds at multiple locations in varied 
patterns near strongly T3 induced genes and this correlates with 
enhancement of T3 response. 

Changes in TRfJ binding events after hormone treatment 

While there were no major changes in overall TR(J distribution 
+/— T3 (Fig. 1), we detected more TRfJ peaks with T3 versus 
untreated cells (6792 versus 5791; 17% increase). We compared 
positions of TR(J binding peaks +/— T3 in more detail. 
Surprisingly, the majority of TRfJ peaks exhibited hormone- 
dependency with only 30% completely unchanged (Fig. S3A). We 
did not observe any specific association of T3-dependent peaks 
and particular functional regions of the genome (Fig. S3B, S3C); 
distributions of ligand dependent binding events reflect overall 
distributions of TRfJ (see Fig. 1). 

To better understand this effect, we inspected hormone- 
dependent TRfJ binding events near target genes. This approach 
revealed only modest changes in TRfJ peak distribution around 
transcription units rather than large scale alterations in TR(J 
binding pattern (see examples in Fig. 4, Fig. S4A). In some 
instances (see LDLR, Fig. 4A and not shown), there were modest 
shifts in TRfJ peak and footprint position. Further, some called 
peaks appeared completely hormone-dependent (see ADSSL1 R2, 
Fig. 4, PDE2A Rl, Fig. S4A and not shown). Such changes were 
not seen at all genes (EPAS1; ACSL5, Fig. S4A). We were able to 
verify 2-fold hormone-dependent changes in TRfJ binding at some 
locations at which TR binding was completely hormone-depen- 
dent by conventional ChIP (Fig. 4B, LDLR; Fig. S4B, PDE2A and 
not shown). Conversely, we found that peaks which appeared 
unchanged after T3 treatment by BioChIP were generally 
unaffected by hormone when assessed by this method (Fig. S4B). 
We therefore suggest that apparent hormone-dependent changes 
in TRfJ binding site distribution (Fig. S3) reflect, at least in part, 
hormone-dependent changes in positions and size of TRfJ binding 
peaks (Discussion). 

TRp Binding Peaks Contain Complex TREs 

Investigation of sequence composition of genomic regions 
bound by TRfJ revealed elements that resembled the typical 
TRE consensus within the peaks. A query of the top 150 bound 
peaks revealed high prevalence of a single TRE half sites 
(TGAGGTCA) (Fig. 5A), distinct from a previous consensus 
derived from analysis of a smaller number of TREs which 
contained two detectable half sites spaced by four bases (Fig. 5B, 
28). Interestingly, this previous study also suggested that the paired 
G residues (TGA GG TCA) were the most important determinants 
of TRfJ binding [30]. While our analysis confirmed that these 
nucleotides are indeed highly represented, it also revealed higher 
representation of other nucleotides within the half site than 
previously denned (TGAGGTCA). 
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Figure 4. Characterization of TFtp binding near induced genes. A. Patterns of TR(! binding depicted at representations of individual target 
gene loci (LDLR, BCL3, NC0R2,ADSSL1 and SOX7). Blue bars represent genomic binding regions, and the vertical red lines represent peaks, as 
classified by QuEST. The horizontal black bars are regions analyzed by ChlP-PCR (locations of primer amplification). Observed binding patterns 
included 5', 3' and intronic binding events, as shown in genomic data tracks (UCSC Genome Browser). Putative regulatory elements, identified 
through sequence analysis of the genomic regions indicated, are depicted below bound regions in which they occur. B. QPCR of ChIP analysis 
confirming DNA binding in regulatory regions of genes. C. Realtime PCR analysis depicting enhancement of transcription of individual loci in Fig. 4A 
by T3 in the presence of TRfi (*P<0.05 by Student's T-Test). 
doi:1 0.1 371 /journal.pone.0081 1 86.g004 



A more directed query, based on the above results and including 
all sequences within TR-bound regions, confirmed the presence of 
canonical half sites in 50% of bound regions (Table 1). A similar 
query with a degenerate half-site sequence (AGGnCA) produced 
positive results for more than 90% of bound regions. Although 
half-sites are common in the genome, there was approximately 2- 
fold enrichment of both consensus and degenerate half-sites in 
TRfi peaks relative to their overall representation. 




Figure 5. Definition of TRE Consensus. A. Consensus sequence 
with similarity to the classic TR|3 binding half-site discovered by analysis 
of top 150 peaks in BioChIP analysis. B. The previously defined TRp 
consensus obtained from analysis of more than 30 published target 
gene regulatory elements is shown for comparison at right. 
doi:1 0.1 371 /journal.pone.0081 1 86.g005 



Close investigation of TRE sequences within TRfi peaks 
revealed that half-sites were usually associated with a second, less 
well conserved, half-site and that organization of these elements 
commonly resembled half-site configurations of classical TREs, 
including DR-4, ER-6, IP-0 and IP-1 (Table 1). Genome wide, 
there was 14-fold enrichment of putative DR-4 elements in TRfi 
binding peaks and 4-5 fold enrichment of ER and IR elements. 
We also detected high representation of canonical DR elements 
with unusual spacing (DR-0 to DR-12) and variations in half-site 
spacing of ERs and IRs, although degeneracy in the 5 ' half site 
can mean multiple interpretations of organization of these TREs 
(not shown). 

To understand whether putative elements identified by com- 
putational approaches were functional TREs, we performed 
further analysis of TRfi binding near the adrenomedullin (adm) 
gene, which is a verified TR target in rodents [58,59]. It is also 
known that adm levels increase in hyperthyroid human patients 
[60,61], but the molecular basis of the latter effect have not been 
previously characterized. We confirmed that T3-dependent 
transcriptional activation of human adm was unabated by the 
translation inhibitor CHX in TRfi-BioChIP cells, proving that it is 
a direct target (Fig. 6 A) and that an adm promoter (— 1 KB) driven 
luciferase reporter also displayed T3 induction after transfection 
into B7B cells (Fig. 6B). There were TRfi binding footprints (blue 
bars) and called peaks (red bars) in the proximal promoter and 
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Table 1. Occurrence of binding motifs in TRR-bound peaks. 





Motif 


Untreated 


Enrichment 


T3-Treated 


Enrichment 


Halfsite 


AGGTCA 


50.32% 


2.10 


48.04% 


2.08 


Halfsite 


AGGnCA 


91.71% 


8.68 


90.75% 


8.51 


DR4 


AGGnCA-4- 
AGGnCA 


1.91% 


17.57 


1.69% 


14.45 


DR4 


AGGnnn-4- 
AGGnCA 


1 1 .36% 


26.00 


1 1 .09% 


25.15 


ER6 


nnnCCT-6- 
AGGnCA 


8.87% 


4.43 


8.01% 


4.26 


IPO 


nnnCCT-0- 
AGGnCA 


3.45% 


1.68 


3.28% 


1.65 


IP1 


nnnCCT-1- 
AGGnCA 


11.81% 


6.10 


10.92% 


5.38 


Negl 


TTTGGG 


43.57% 


1.76 


42.37% 


1.75 


Neg2 


CCCCTCAGGCGC 


0.02% 


2.09 


0.01% 


1.87 





doi:1 0.1 371 /journal.pone.0081 1 86.t001 



immediately downstream of the gene (Fig. 6C) and we verified 
these binding events with conventional ChIP using indicated 
primers (R1-R4; Fig. 6D). Use of a custom position-weighted 
matrix-based comparison with previously characterized TR 
regulatory sites revealed two probable TREs within R2, a DR-4 
element and a DR-6 element that lay within 1 KB of the TSS 
(Fig. 6C). Binding of both elements to TRP was confirmed by 
EMSA (Fig. 6E) and both robustiy activated transcription from a 
luciferase reporter (Fig. 6F). Thus, T3 induction of human adm 
involves TRP binding to TREs in the proximal promoter, along 
with other sites. We were routinely able to identify TREs in other 
TR(3 peaks and verify function (see Fig. 7, Figs. 4 and SI and 
manuscripts in preparation). We conclude that TRP is commonly 
recruited to DNA through interactions with TREs that contain 
canonical or variant AGGTCA half sites, albeit with diversity in 
spacing and orientation. 

Finally, we considered representation of putative negative 
response elements, including "TTTGGG" and 
"CCCCTCAGGCGC" [35,36,39] in TRP binding sites 
(Table 1). Both elements were enriched in TRP peaks relative to 
their overall genomic representation (1.7 and 1.9 fold respectively). 
Investigation of TRP peaks that were proximal to a highly 
negatively regulated gene (PDE2A see Fig. 3), however, failed to 
reveal possible nTREs, although we did observe classical TRE-like 
sequence elements that resembled DR-4 sites (Table SI). Further, 
consideration of peaks detected near other negatively regulated 
genes (SOX9 etc.) that not meet cutoffs indicated in Figure 3 failed 
to reveal any nTRE-like sequences in most cases, the sole 
exception being C90RF169 (Table S6). Thus, hormone-depen- 
dent suppression of TR target genes in HepG2 cells is not 
commonly associated with detectable TRP binding to nearby 
nTREs (Discussion). 

Intergenic TRP Binding Sites are functional TREs 

Since more than 30% of TRP binding events occurred in 
intergenic regions (Fig. 1A), similar to other NRs, we tested the 
ability of some of these sites to regulate transcription. We verified 
TRP binding to three intergenic peaks with ChIP (Fig. 7A, 7B). 
We also used ChIP to confirm that T3 induced an activating 
histone acetylation near each peak (H3Ac, Fig. 7B). We identified 
putative TREs within each peak (Fig. 7A) and confirmed activity 
by EMSA (not shown) and assays of luciferase activity (Fig. 7C). 



Thus, intergenic binding events can induce an active chromatin 
state in response to T3 and co-localize with functional TREs. 

Colocalization of TR and TF Binding Sites 

Finally, we determined whether TRP might associate with 
heterologous TFs in composite modules. We examined TRP peaks 
for over-representation of TF binding sites, using a comprehensive 
query of regions with known TF consensus sequences in the 
JASPAR database [54], with Seqpos (http://cistrome.org) and 
these are listed in Table 2, with the most commonly associated 
factors listed at top. 

We identified a variety of NR binding sites, including RAR, 
VDR, PPAR, RORot, ERoc, ERP and Nurrl, all of which consist 
of canonical AGGTCA half-sites of distinct spacing (e.g., RAR 
binds to DR-5 elements, VDR binds DR-3 and PPARs bind DR- 
1). While we cannot exclude the possibility that TRs associate with 
these NRs, the simplest interpretation of this finding is that these 
elements are atypical TREs, consistent with promiscuity of TRP 
DNA binding, see Discussion. 

We also observed co-localization of TRP binding sites with 
consensus binding sites for several heterologous TFs. These 
include AP-1 (jun), CREB, p53 and CTCF, all of which have 
been previously shown to associate with TRP (Table 2). Addi- 
tionally, we detected association of TRP with NFE2L2 binding 
sites (also known as NRF2), involved in regulation of the 
antioxidant response pathway, AP2 and others. There were no 
major changes in preference of TRE/TF binding site co- 
localization near TRP peaks +/— T3 (Fig. S2). Thus, our findings 
support previous studies which suggest that TR cooperates with a 
small subset of heterologous TFs that includes AP-1 and CTCF 
and also identify new candidates for consideration. 

Discussion 

In this study, we have assessed genomic distributions of TRP in 
a human liver cell line to understand how TRP binding patterns 
are related to T3 and TR-dependent changes in gene expression. 
We used HepG2 cells that express exogenous BLRP-tagged TRP 
and enzymatic machinery required for biotinylation (BirA) to 
examine TRP binding [44] because we and others have found that 
HepG2 cells display highly reproducible and widespread T3- 
responses at physiologically relevant genes [54,62]. Stable TR 
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Figure 6. Links between TR Binding and Adm Transcription. A. Graph showing results of realtime PCR analysis of adm transcription in the B7B 
cells after six hours of T3 treatment +/— 10 |ig/ml CHX cotreatment of B7B cells. B. Patterns of TRp binding peaks at the adm locus (UCSC Genome 
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Browser), in similar format to Fig. 4. TRp binding events clustered into four regions (R1, R2, R3, R4), upstream and downstream of this transcript, as 
well as a substantial amount of binding immediately proximal to the transcriptional start site. C. Binding of TRp was confirmed by realtime ChIP PCR 
analysis in B7B cells at the regions indicated (ChIP primers are depicted by horizontal bars in B). D. The proximal promoter region of adm 
(corresponding to R2) conferred T3-dependent increases in luciferase activity upon a standard reporter after transfection into B7B. E. Results of gel 
shift confirming direct TRp binding to two putative response elements, designated TRE-1 and TRE-2 that were found in R2 at positions marked in 
Fig. 6B. Individual lanes show shifts obtained with elements and RXRoc-TRp +/- competitor DNA or mutated versions of both elements. F. Luciferase 
reporter assays confirming that TRE-1 and TRE-2 confer T3 responsiveness on a reporter gene. (*P<0.05 by Student's T-Test). 
doi:1 0.1 371 /journal.pone.0081 1 86.g006 



expression is necessary because TR mRNA levels are greatly 
reduced in liver cell lines and primary cultures versus native [57]. 
Given this requirement, we selected the BioChiP system because it 
allows straightforward, high stringency purification of TR(3 [44]. 
Interestingly, many of our findings resemble those obtained with 



exogenous TRs in a mouse neural cell line suggesting that 
combined results of these approaches are revealing at least some 
fundamental principles of TR genomic distributions and actions 
[7,19]. 
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Figure 7. Intergenic binding events. A. Three intergenic binding peaks were selected and analyzed for the presence of recognizable TR binding 
motifs (sequences of motifs listed). B. Results of qPCR ChIP analysis confirming binding of TRp to the intergenic regions depicted in Fig. 7A (top 
panels) and induction of H3 acetylation near sites (bottom panels). C. Results of luciferase reporter assays, with indicated constructs containing 
intergenic elements described in Fig. 7A, confirming that each element confers T3 induction. (*P<0.05 by Student's T-Test). 
doi:1 0.1 371 /journal.pone.0081 1 86.g007 
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Table 2. Occurrence of previously identified binding motifs in 
TRfl-bound peaks. 



ID 


Name 


Family 


Z-Score 


MA0099 


JUN 


BZIP 


-32.85 


MA0160 


NURR1 


NHR 


-27.44 


MA0018 


CREB1 


BZIP 


-26.47 


MA0071 


RORA 


NHR 


-27.6 


MA0159 


RARA 


NHR 


-15.27 


MA0112 


ESRA 


NHR 


-15.21 


MA0150 


NFE2L2 


BZIP 


-13.43 


MA0258 


ESRB 


NHR 


-13.3 


MA0115 


LXRB 


NHR 


-12.17 


MA0043 


HLF 


BZIP 


- 1 0.38 


MA0003 


TFAP2A 


AP2 


-10.37 


MA0074 


VDR 


NHR 


-9.72 


MA0017 


COUPTF1 


NHR 


-8.62 


MA00139 


CTCF 


GO 


-7.91 


MA0066 


PPARG 


NHR 


-7.76 
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Our findings indicate that mechanisms of T3-dependent gene 
induction often conform to prevailing models of TRP action. 
While TRp 1 distribution resembles that of other NRs, with TRP 
peaks spread widely through the genome and often far from 
obvious TR target genes, we nevertheless detect striking enrich- 
ment of TRP at proximal promoters, the TSS and immediate 
downstream regions of transcribed genes (Figs. IB, C) and more 
than 45% of T3 inducible genes are associated with TRP peaks 
(Fig. 2B). Further, >90% of TRP peaks contain a TRE half-site, 
there is significant enrichment of DR-4 elements and other known 
TRE configurations and we commonly identify putative TREs 
within bound regions (see Figs. 4, 6 and 7 and Table S1,S6) and 
are able verify TRE function (Figs. 5 and 7 and manuscripts in 
preparation). 

More surprisingly, we failed to detect obvious TRP enrichment 
near T3 repressed genes. Further, while we do observe enrichment 
of putative nTREs in TRf} peaks, we have been unable to link 
specific nTREs to nearby highly repressed targets. Both findings 
appear surprising in light of previous results which show that TRs 
are directly recruited to negatively regulated genes such as [34— 
36], but we think that our results do not contradict these findings 
and instead point to the existence of alternate mechanisms of TR 
gene repression that are active in HepG2 cells. We previously 
found that large proportions of negatively regulated genes are 
inhibited by CHX in HepG2 cells, implying that new protein 
synthesis is often required for transcriptional repression in this 
context, and we therefore suspect that many individual instances 
of T3 repression in HepG2 cells could involve secondary 
responses. Other phenomena may also be at play, including 
TR/T3-dependent remodeling of coregulator complexes that 
interact with multiple genes. Clearly, this issue will require further 
investigation and it is not clear whether a similar preponderance of 
indirect effects will also be a feature of T3-dependent negative 
regulation in vivo or whether this is restricted to the HepG2 system. 
We also failed to detect significant TR(3 enrichment near genes 
that respond to unliganded TRs. We think that this could also 
reflect a preponderance of indirect response to active unliganded 
TRs. 



While TR actions at many T3 induced genes often appear to be 
mediated by proximal binding events, our findings emphasize that 
mechanisms of T3 induction are not always straightforward; 
around 50% of T3 induced genes lack proximal TRf3 binding 
peaks, raising the possibility that such instances are explained by 
either long-range regulatory events or local binding events which 
are too ephemeral for detection through the currently employed 
method. While we expect that the latter possibility will be true in 
some cases, we surveyed several proposed TRP binding sites that 
were not detected in our ChlPseq studies, including within the 
pckl promoter, and have concluded that functional TREs are not 
present and that we did not overlook TRP binding in these cases 
(not shown). Another possibility, mentioned for T3 repression, is 
that some T3 induced genes are indirect TR targets that respond 
to primary TR-dependent changes in levels of other TFs. We 
suspect that this only applies to a small number of examples of T3 
induced genes that lack TRP binding sites; genome wide 
assessment of CHX sensitivity in HepG2 indicates that most T3 
induced genes are direct TR targets [54]. We therefore instead 
favor the idea that this subset of T3 induced genes is subject to 
extremely far distant regulation by functional intergenic TRP 
binding sites [4, 1 7] . While it is intriguing to suggest that putative 
intergenic TREs, similar to those in Fig. 7, are responsible for 
these long range effects, this idea must be treated with caution. 
The fact that these elements are not near known genes but activate 
luciferase from a proximal location suggests that they have the 
capacity to be functional but their true role in the context of the 
whole genome is unclear. While the increase in acetylation seen 
nearby these locations with T3 is interesting, it is also conceivable 
that TRP binds adventitiously to elements that happen to lie 
within "open chromatin" rather than to true functional elements 
in this system. 

The ChlPseq approach confirms and extends previous sugges- 
tions that TRE organization is variable [63,64]. We have been 
able to define a consensus TRE half-site that is commonly 
represented with TRP peaks, found that consensus half-sites 
(AGGnCA) occur in >90% of TRP peaks and observed 
enrichment of classical DR-4 elements and, to a lesser extent, 
other canonical TREs. However, total representation of DR-4, IR- 
0 and ER-6 elements is much lower than total TRE half-site 
representation (Table 1). This discrepancy is a consequence of 
high representation of TRE-like sequences with non-canonical 
half-site spacing within TRP peaks (Table 2) and probably also 
explains why we detected large numbers of binding sites for other 
NRs that interact with AGGTCA half-sites. This raises obvious 
questions about: i) the structural basis of flexibility of TR/TRE 
recognition, ii) connections between promiscuous DNA recogni- 
tion and gene-specific variations in TR activity and iii) possibilities 
for NR/TR cross-talk in gene regulation. Our study does reveal 
differences with the canonical TRE site identified by Chatonnet 
and coworkers, who observed a DR-4 like consensus in neural cells 
rather than the half-site seen here [19]. We do not understand 
whether this difference is related to technical aspects of our study 
versus that of Chatennet, including TR expression levels, or 
whether there are true differences in response element recognition 
patterns in different cell types. This issue will also require further 
investigation. 

Our findings also reveal unexpected features of TR binding site 
architecture. T3-induced genes are often associated with clusters of 
TRP peaks rather than single elements, and there is enrichment of 
TRP binding both 5' and 3' of transcribed genes and binding sites 
can be seen within untranscribed regions and within introns. It is 
not clear whether TRs play distinct roles in gene expression when 
bound to different locations with respect to the transcription unit; 
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this issue will require further investigation. It is also intriguing that 
computerized analysis of consensus TREs throughout the genome 
indicates that they are not obviously over-represented near TR 
target genes, implying that actual TR(5 binding events are 
dependent on other factors, possibly local states of chromatin 
modification or binding events of partner TFs. It will be important 
to understand mechanisms that underlie TR(5 binding site 
selection within the milieu of living cells. 

A large proportion of TR(5 peaks exhibit some apparent T3- 
dependency, but we do not believe that this reflects large scale 
redistribution of TRs in response to hormone. Overall genomic 
localization of TRs is similar in the absence and presence of T3 
(Fig. S3) and reanalysis of data with altered stringency of peak 
calling did not change our conclusion that overall TR distribution 
does not change after hormone treatment (not shown). Close 
investigation of TRP peaks near target genes revealed relatively 
modest alterations in TR(3 peak position and footprint size rather 
than large scale appearance or disappearance of TRs from the 
vicinity of target genes (Fig. 4 and Fig. S4). We found some 
examples in which apparent changes in footprint size could not be 
verified with standard ChIP (see adm footprints in Fig. 5). More 
commonly, however, we were able to verify at least some degree of 
hormone-dependency of TR(3 binding at selected peaks with 
conventional ChIP. Thus, we favor the idea that there are modest 
redistributions in TRP binding after hormone treatment and think 
that this effect accounts for the large apparent change in TR(3 
peak distribution after hormone treatment. 

Functional significance of verifiable hormone-dependent chang- 
es, if any, remains unclear. T3 may promote relocalization of TRP 
from inactive DNA pools near a negatively regulated target gene 
to nearby functional regulatory elements [65]. Our observations 
suggest, however, that T3-dependent changes in TR(3 footprint 
and peak position are complex and gene-specific with no obvious 
pattern. We recently showed that hormone-dependent TRP 
binding to a TRE within the glucose-6-phosphatase promoter 
requires TR(5 interactions with a gene specific cofactor, the 
NAD+-dependent deacetylase Sirtuin 1 (SIRT1 [66]). Clearly, 
more investigation will be needed to understand the potential 
importance of this unexpected phenomenon in the context of T3- 
dependent TR(3 binding. 

Finally, our results allow us to define possibilities for TR cross- 
talk with other TFs. Fewer than 10% of TRP peaks lack 
identifiable TREs, suggesting that DNA independent TRP 
recruitment that exclusively involves contacts with heterologous 
TFs such as AP-1 is relatively rare. However, we do detect strong 
enrichment of binding sites for TFs that are known to interact with 
TRs, including AP-1 [37,38], CTCF [39] and p53 [40], and for 
other TFs that have not previously been proposed to cooperate 
with TRs, including NRF2 and AP2. It will be interesting to see 
whether TR interacts with these factors in composite modules and 
examine roles of complex elements in T3 response and cross-talk 
of TRs with other signaling pathways. 

In summary, our studies indicate that T3 gene regulation 
commonly involves proximal TRP binding events near target 
genes, but also reveals striking variability in TRE position, distance 
from the target gene, TRE half-site organization, sequence and 
spacing and interactions with heterologous TFs. We and others 
have previously described large variations in precise mechanisms 
of TR action at different target genes [54,57,67]. It is interesting to 
speculate that these differential effects are linked to variations in 
genomic context of TR(3 binding sites described here. 



Supporting Information 

Figure SI BioChlP-modified TR0 constructs exhibit 
similar activity to endogenous TRp. A. BLRP-tagged TRP 
was expressed in HepG2 cells at a level 8-10 fold higher than that 
endogenously expressed in parental HepG2 cells, as assessed by 
realtime PGR, normalized to 18S RNA. B. TR(3 protein levels in 
B7B cells (lanes 5-8, 20 u.g cellular protein per lane) were found to 
be substantially less than that expressed in mouse liver (lane 9, 
20 |xg protein per lane) as assessed by western blot. C. Modes of 
DNA binding were assessed by EMSA assays, as indicated. B7B 
cells elicited mobility shift of a DR4 element (left panel, lanes 6-9), 
similar to that of previously characterized [54] Flag-tagged TRP 
cell lysates ("HG2 FlagB", lanes 10-13). Supershift of complexes 
with RXRa or TRP antibodies indicated the predominandy 
heterodimeric composition of complexes (left panel, upper band). 
Analysis of B7B cell lysates, in comparison to in vitro-expressed 
RXRa and TRP protein alone (middle panel, lanes 2&3) or 
combined lysates (lane 4) confirmed a predominandy heterodi- 
meric mode of binding to DR4 elements (lanes 5-7), similar to 
Flag TRP lysates (lanes 8-10). Analysis of T3-treated cell lysates 
revealed a characteristic change in mobility, as compared with 
untreated lysates, similar in B7B (right panel, lanes 5-7) and HG2 
FlagB lysates (lanes 8-10). D. Transcriptional activation of DR4 
and ER6 elements were quantified in comparison to controls in 
luciferase reporter assays in HepG2 cells (left panel) or TRP- 
expressing HepG2 cells (right panel). E. Expression of defined 
TRP target genes was quantified by realtime PCR of target genes 
in indicated cell models, normalized to 18S RNA. F. Protein 
expression of the TRP-regulated gene LDLR was assessed by 
western blot in TRP-expressing HepG2 cells, showing substan- 
tially increased protein levels after T3 activation of TRp. 
(TIF) 

Figure S2 Comparison of B7B cells and HepG2-TRp 
(Flag) cells. A. Expression of BLRP-tagged TRP or Flag-tagged 
TRP was 8-10 fold higher than endogenous TRP expression in 
HepG2 parental cells and was not affected by treatment with 
100 nM T3 for 8 hours, as assessed by realtime PCR, normalized 
to 18S RNA.Expression of TRa was similar in all samples. B. 
Expression of defined TRP target genes was induced with 8 hour 
T3 treatment to a similar extent in BLRP-tagged TRP and Flag- 
tagged TRP cells, as assessed by realtime PCR of indicated targets, 
normalized to 18S RNA. C. A DR4 luciferase reporter construct 
showed a similar transcriptional activation profile after treatment 
of BLRP-tagged TRP and Flag-tagged TRP cells; little activation 
was observed in HepG2 parental cells. 
(TIF) 

Figure S3 Summary of binding events in TrP-expressing 
cells. A. TRP binding peaks were mapped to genomic regions, 
using the methods described. Approximately 30% of binding to 
sites occurred only in the presence or absence of T3, as depicted in 
the Venn diagram shown. B. Percentages of total binding events 
within specific regions, as assessed with the CEAS analysis 
program [68,69], are shown, including those occurring exclusively 
in untreated samples, or exclusively in T3-treated samples are 
indicated. C. Percentages of binding events, which were 
exclusively observed in untreated or T3-treated samples, was 
mapped to the genomic regions indicated, using the CEAS 
analysis program, and depicted in the charts as labeled. 
(TIF) 

Figure S4 Characterization of TRp binding near in- 
duced genes. A. Patterns of TRP binding were depicted at 
individual target gene loci, represented using the same format as 
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Fig. 4. Observed binding patterns included 5', 3' and intronic 
binding events, as shown in genomic data tracks (UCSC Genome 
Browser). B. Patterns of TR(3 binding at the regions indicated with 
ChIP, analyzed by QPCR. (*P<0.05 by Student's T-Test). 
(TIF) 

Table SI Response elements characterized in reported 
data. Hormone response elements are listed by their associated 
genes (Column 1). The DNA sequence of each element is 
provided, with spaces to indicate half-site position (Column 2), 
and the assays used to evaluate the elements described (Column 3). 
(PDF) 

Table S2 Primers used for chromatin immunoprecipi- 
tation. TRB binding at indicated peaks was assessed by realtime 
PGR of chromatin IP-enriched samples with indicated primer sets. 
Primer sets were designed to specifically amplify indicated 
genomic regions. 
(DOC) 

Table S3 Perl script used for analysis of bound 
sequences. Sequence composition of TRB-binding genomic 
regions was analyzed with the included Perl script, confirmed to 
run on Ubuntu Standard Distribution 11.1, in addition to the 
analysis procedures described in text. This script analyzes two 
BED-formatted input files by evaluating the length-normalized 
bidirectional site occurrences. 
(PDF) 

Table S4 Primers used for realtime PCR-based analysis 
of gene expression. Gene expression was assessed by realtime 
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